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We study the drag on an object moving through a trapped 
Bose-Einstein condensate, and show that finite compressibil- 
ity leads to a mutual drag, which is subsequently suppressed 
by the formation of a vortex ring. 

PACS numbers: 03.75.Fi, 67.40.Vs, 67.57.De 

The appearance of vortices in quantum fluids has im- 
portant consequences for many physical systems |IJ|^] . In 
particular, there appears to be a direct link between vor- 
tex formation and the breakdown of superfluidity in liq- 
uid helium However, a quantative comparision be- 
tween theory and experiment is impeded by the com- 
plexity of this many-body system. 

The recent observation of Bose-Einstein Condensation 
(BEC) in trapped alkali gases 0], on the other hand, 
has yielded an ideal testing ground for many-body theo- 
ries. Dilute condensates can now be produced at temper- 
atures far below the BEC transition, allowing an accu- 
rate description by the Gross-Pitaevskii (GP) equation. 
By studying this equation, insight can be gained into the 
relationship between vortex formation and drag in super- 
fluid systems. For homogeneous fluid flow past an obsta- 
cle, numerical studies show that below a critical velocity 
the fluid exerts no net force on the object ||, whereas 
above this velocity vortices are nucleated, leading to a 
pressure imbalance which creates drag Q . 

The picture is less clear, however, in trapped Bose 
gases, which are inhomogeneous and of finite size. Simu- 
lations have shown that dragging an object, created by a 
far-detuned laser beam, through an inhomogeneous con- 
densate produces compression waves and a series of vor- 
tex pairs 0. The formation of these structures will in- 
fluence the dynamics of the condensates and the drag on 
the object. An attractive system for studying drag is the 
multi-component condensate ^| , where the magnetic sen- 
sitivity of the different hyperfine states allows the fluid 
components to be intermingled or spatially separated. In 
this paper, we simulate actively dragging a small conden- 
sate, the 'object', through a larger, less tightly confined 
condensate, the 'fluid'. We find that a mutual drag arises 
due to deformation of the compressible object, which is 
relieved by the formation of vortex rings in the fluid. The 
time dependence of the mutual drag provides a signature 
of ring formation. This process is analogous to ring nu- 
cleation by moving ions in superfluid 4 He ||, albeit on a 
different scale. 

The numerical methods follow closely those employed 
in our previous work ]7|,^0|. The coupled GP equations 



for the condensate wavefunctions, ipi(r,t), in the x — y 
plane can be written as: 

idtik = [-V 2 +Vi + C (c^| 2 + IV^I 2 ) - Mi] V'i, 

(1) 

where i, j — I, 2 (i ^ j), ^ are the chemical potentials at 
equilibrium, Vi are the trap potentials, on are the ratios of 
the scattering lengths, and C is the nonlinear coefficient 
p| . We have performed 2D, cylindrically-symmetric 3D, 
and full Cartesian 3D simulations ]l2| . 

To demonstrate that coupled GP equations give an ac- 
curate description of the formation of structure in the 
condensate, we begin by modelling an experiment on 
component separation p3| . We assume that a conden- 
sate initially in state |1) with trapping potential V\ = 
(x 2 + e 2 y 2 ) /4, is subject to an interaction which transfers 
50% of the population into state |2), which experiences 
a displaced potential Vi — (x 2 + e 2 (y — z/o) 2 )/4, where 
?/o is the offset Q. The displaced potential causes the 
two components to separate and subsequently oscillate 
at large amplitudes. Nonlinear mixing |l5| ] leads to ex- 
citation of high-order modes, which appear as structure 
in the density profiles |l3j. Fig. 1(a) displays a density 
cross-section at t = 8.52 (corresponding to 65 ms). How- 
ever, in the experiment, most of the structure had dis- 
appeared by this time. This disparity may be due to 
Landau damping at finite temperatures, where conden- 
sate excitations are absorbed by the thermal cloud. To 
simulate this effect, we propagate the GP equations (|l|) 
in complex time t — ► (1 + iA)i p6| . Components with 
an energy, E, decay exponentially at a rate proportional 
to the damping constant, A < ]l7| ], and to E — 
in agreement with theory for the collisionless low-energy 
regime Jl8{ . Higher frequency excitations are preferen- 
tially damped, reducing the density variation within the 
condensate (Fig. 1(b)). 

Fig. 2 shows the center-of-mass position of the two 
condensates (Yi — JJ y\tpi\ 2 dxdy) as a function of time. 
Population of high-energy excitations results in damp- 
ing of the dipole mode, in agreement with previous work 
fll9[ |. However, this mechanism alone is insufficient to ex- 
plain the observed experimental damping. Propagation 
in complex time results in further damping, but at a rela- 
tively slow rate. We also investigate the effect of an offset 
in the a;-direction, which may arise as a consequence of 
the experimental geometry |20j. In this case, mixing be- 
tween the fast (y) and slow (x) vibrational modes |2l| ] 
could be more significant than finite temperature effects 
in explaining the observed damping (see Fig. 2), though 
it is likely that a combination of these mechanisms is 
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responsible. 
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FIG. 1. Cross-sections through 2D density profiles, |^i| 2 
(dashed) and |t/>2| 2 (solid), at x — and t = 8.52, showing 
results of (a) undamped, and (b) damped (A = —0.025) prop- 
agation. 
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FIG. 2. Center-of-mass oscillations, Yi(t), in |1) (top) and 
|2), for undamped (solid) and damped A = —0.025 (dashed) 
propagation. The dot-dashed lines show the effect on un- 
damped propagation of an additional offset, xo = 0.240, 
where decay of the oscillations are more marked than even 
a large Landau damping. 

As is apparent in Fig. 2, over short timescales (t < 
2) thermal damping is negligible and the undamped GP 
equations provide a reliable model. We now turn to the 
situation where one condensate flows through the other, 
in direct analogy with an object moving through a fluid. 
This may be realized experimentally by using a magnetic 
trap to confine atoms in state |2), whilst an optical dipole 
trap, moving with relative velocity —v, loosely confines 
atoms in a magnetically insensitive level |1). We employ 
the coupled equations (Q), with V\ — (x 2 + y 2 )/4 and 
V 2 = V\ + {x 2 + (y - vt) 2 ), and on = 1. The initial 
state is found using imaginary time propagation |2^j . The 
repulsive mean field arising from the 'object' (|2)) creates 



a local minimum in the center of the density profile of the 
'fluid' (|1)) (^3|. The depth of the minimum depends on 
the interaction strength C, and on the fraction of atoms 
in the fluid, /. 
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FIG. 3. Cross-sections through 2D density profiles at 
x = 0. Condensate \2) (the 'object', solid line) moves through 
|1) (the 'background fluid', dashed), due to displacement 
of trap potentials at constant velocity, v — 3; C = 1050, 
/ = 0.95 (where / is the fraction of atoms in |1)). Two 
time-frames are shown: (a) t — 0.75, where drag arises from 
a process analogous to phonon emission from the accelerat- 
ing object, (b) t — 1.3, where deformation of the object and 
surrounding fluid leads to an additional drag. 

Displacement of the object potential (V2) at t > in- 
duces motion of the object, which leads to the minimum 
in the background fluid becoming progressively deeper at 
a rate which increases with v (see Fig. 3). When the den- 
sity minimum reaches zero, it evolves into a vortex ring. 
In addition, motion of the object creates a response in 
the fluid, implying a back-action, or drag, on the object. 
The drag may be studied through the center-of-mass ac- 
celeration, given by: 



Ri(t) = -2 J \^(r,t)\ 2 Vfi t (r,t) d' 



r. 



(2) 



For steady flow of an homogeneous condensate without 
vorticity, the integral vanishes, implying superfluidity 
p4j . However, for a finite inhomogeneous system there 
is no steady state and a different criterion for superfluid- 
ity should be sought. There are four energy terms that 
contribute to the acceleration; however, the kinetic and 
self-interaction terms are found to be negligible, leaving 
only: 



and, 



oy 



d 



y,,nnAI) : -2C \^\ 2 — \^\ 2 dxdy, 



(3) 



(4) 
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which correspond to the contributions from the trap po- 
tentials and the interaction between condensates, respec- 
tively. 

The mutual acceleration (^) is particularly interesting 
in addressing the issue of drag. Fig. 4 shows the calcu- 
lated mutual acceleration as a function of time. At later 
times (t > 2) the condensates repel due to the inhomo- 
geneity in the background fluid. For t < 2 the accel- 
eration is negative, implying an effective attraction, or 
equivalently a drag on the object. 



system H|). However, due to the finite size of the back- 
ground fluid there is no steady flow condition, and drag 
is produced even at low velocities. 
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FIG. 4. Plot of y a ,aut(t) for C = 600, / = 0.83; 
v = 1,2,3,4. The 'drag' is evaluated from the minima of 
the curves, which increases in magnitude as the velocity rises. 
Inset shows a similar plot (C = 1050, / = 0.95) for v = 0.25 
to v = 1.75 in steps of 0.25. At low v the local minima arise 
due to the finite response of the background fluid, |1). At 
higher v an additional minimum appears at t > 1 due to the 
compressibility of the object. The position of each minimum 
coincides with the moment of vortex formation. 

For t < 1 and low velocities, Fig. 4 (inset), the force 
results from the slow response of the fluid to the object 
acceleration. The object moves to the front of the poten- 
tial well it creates in the mean field of the fluid, resulting 
in a restoring force which persists until the fluid can re- 
spond (Fig. 3(a)). The maximum attractive drag (i.e. 
the minima in Fig. 4) is plotted as a function of velocity 
in Fig. 5. This process, which is equivalent to phonon 
emission by an accelerating object, is responsible for the 
linear section of the drag curve. 

For velocities near the speed of sound in the object, 
c 2 = ^/2C|V'2| 2 , the object begins to deform, with the 
result that the overlap between the two fluids is reduced 
behind and enhanced in front (Fig. 3(b)). This increases 
the drag, producing the additional minima at t > 1 in 
Fig. 4 (inset). Onset of this process is rapid with in- 
creasing velocity, giving rise to a sharp transition in the 
drag curve. This transition can be compared to that 
between superflow and normal flow in the homogeneous 
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FIG. 5. Peak mutual attraction (drag) as a function of 
velocity. Curves are shown for: * C = 600, / = 0.83; A 
C = 1100, / = 0.91; □ C = 1200, / = 0.83; • C = 1050, 
/ = 0.95. The curves show linear dependence of drag at 
low velocities, and enhanced drag at high velocity due to the 
compressibility of the condensates. The magnitude of the 
additional drag is enhanced for highly compressible object 
condensates i.e. when 1 — / is small, or when the mean-field 
interaction between condensates is large. 

The drag force increases with time until the local min- 
imum in the background fluid reaches zero (Fig. 3(b)), 
from which a vortex ring is formed. Subsequent expan- 
sion of the ring results in the condensate minimum being 
filled, thereby decreasing the pressure imbalance across 
the object. The object returns to a more symmetric 
shape, and the drag decays, as is apparent in Fig. 4. 
Hence, vortex formation in this context tends to limit or 
reduce the drag. 

In summary, we simulate a scheme where one conden- 
sate is pulled through another, and identify mutual drag 
effects, and vortex ring creation which acts to suppress 
drag. Although vortex rings would be difficult to detect 
directly, we have shown that they produce a significant 
change in the mutual drag which could be observable. 
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